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We introduce a procedure to generate scattering states which display trajectory-like wave func- 
tion patterns in wave transport through complex scatterers. These deterministic scattering states 
feature the dual property of being eigenstates to the Wigner-Smith time-delay matrix Q and to 
the transmission matrix t't with classical (noiseless) transmission eigenvalues close to or 1. Our 
procedure to create such beam-like states is based solely on the scattering matrix and successfully 
tested numerically for regular, chaotic and disordered cavities. These results pave the way for the 
experimental realization of highly collimated wave fronts in transport through complex media with 
possible applications like secure and low-power communication. 
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The scattering of waves through complex systems is a 
central subject in physics occurring on a variety of length 
and time scales. Coherent electron transport through 
mesoscopic systems, light transmission through optical 
devices as well as all matters related to room acoustics 
are just a few examples of this kind. Recently, enor- 
mous experimental progress has been made in the abil- 
ity to determine the system-specific scattering matrix of 
such complex systems either explicitly [l[ or implicitly 
by methods like adaptive wave-front shaping Q and op- 
tical phase conjugation These advances have led to 
spectacular results for complex scatterers which could be 
made transparent [H, El or put to use for focusing an in- 
cident wave on a spot size below the diffraction limit 0] ■ 

Common to all such applications is the aim to employ 
the information stored in the scattering matrix to create 
scattering states with specific properties. A very funda- 
mental property a scattering state can have is to follow 
the particle- like bouncing pattern of a classical trajectory 
throughout the entire scattering process Such "clas- 
sical" scattering states play a key role for the wave-to- 
particle crossover, for the emergence of geometrical optics 
out of physical optics and for the breakdown of univer- 
sality in coherent transport [6MToj] . However interesting 
these classical states may be, in complex scattering ge- 
ometries they turn out to be as elusive as the proverbial 
needle in a haystack. 

In this Letter we propose an operational procedure to 
generate such states explicitly. Our approach is illus- 
trated with the example of a two-dimensional rectangular 
cavity through which waves can be scattered by two leads 
attached to the left and right (see Fig.Q}. With each lead 
carrying TV open modes the (2N x 27V)-dimensional uni- 
tary scattering matrix of this device has the form, 



S 



t' 



(1) 



where each of the four blocks contains N x N complex 
elements for the energy dependent transmission (t) and 
reflection (r) amplitudes [(unprimed) primed amplitudes 



designate injection from the (left) right lead]. The to- 
tal transmission T through this resonant cavity is given 
as T = Tr(th) = £„ =1 r„, where the r„ e [0, 1] arc the 
real transmission eigenvalues of the hermitian matrix t't. 
Among the associated eigenstates |r) those with eigenval- 
ues close to r = or r = 1 are termed "noiseless states" as 
they feature a vanishing contribution to electronic shot 
noise Since all of the desired "classical" states 

with a trajectory-like bouncing pattern must have such 
deterministic values of transmission they are all part of 
a highly degenerate noiseless subspace associated with 
t = 0,1. Consider in Fig. QJi,b two randomly chosen 
states with r > 0.99 from this subspace, calculated with 
the Modular Recursive Green's Function Method 
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FIG. 1: (Color online) Scattering through a rectangular cav- 
ity (flux injected through the left lead of width d): (a),(b) 
Wave function densities of transmission eigenstates |r) of t t 
with similar transmission eigenvalues r > 0.99 but different 
wave numbers: (a) k = 5.5n/d, (b) k — 75.5n/d. (c) Classical 
surface of section, recorded for trajectories which enter at the 
left lead mouth with vertical position y and transverse mo- 
mentum p„. The largest of the transmission/reflection bands 
(black/white) are labeled by T1,T2/R1,R2 (bands which are 
equivalent in an extended zone scheme are given the same 
label), (d) Husimi distributions of the states shown in (a) 
(upper panel) and in (b) (lower panel). The size (area) of 
the Planck cell h is indicated by dashed black frames and the 
underlying classical phase space is shown in gray. 
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The first such state (see Fig. QJ,) was calculated at 
a low wave number where only N — 5 lead modes are 
open. To understand the composition of classical trajec- 
tories contributing to this state we evaluate the Poincare 
surface of section (PSS) at the entrance lead junction 
(see Fig. QJ). With transmitted (reflected) trajectories 
being shown in black (white) the PSS features a banded 
pattern with the individual bands being made up of bun- 
dels of trajectories which all have equivalent bouncing 
patterns [12J. The contributions of different phase space 
bands to the state |r) in Fig. [TJ, are revealed by compar- 
ing its quantum phase space distribution (Husimi func- 
tion) H(y,p y ) = \(T\y,p y )\ 2 with the PSS, where \y,p y ) 
is a minimum uncertainty state at y,p y - The correspond- 
ing plot in Fig. QJi (upper panel) shows that the banded 
structure of the PSS is not resolved by this state — in line 
with the fact that the individual areas of the largest phase 
space bands are all smaller than the Planck constant h, 
i.e., the lower resolution limit in wave scattering. Rather, 
this state is composed of many interfering contributions 
from both transmitting and reflecting bands. 

The situation is different when, for smaller wave 
lengths, the size of the Planck cell is well below the size 
of the largest phase space bands Consider, e.g., 

the noiseless transmission eigenstate |r) with N = 75 
shown in Fig. [TJd. We find that the Husimi distribution 
of this state [see Fig. [TJi (lower panel)] is entirely located 
on transmission bands, indicating that full transmission 
is reached here by resolving the classical phase space. 
Since, however, more than one transmission band con- 
tribute (mostly Tl and T2) clear signatures of classical 
bouncing patterns are still absent in the corresponding 
scattering wave function (Fig.QJi). This is a result of the 
indiscriminate mixing of states in the degenerate noise- 
less subspace. This problem may be circumvented by 
explicitly constructing scattering states which lie on in- 
dividual phase space bands |9| . In a real experiment such 
a protocol, however, meets the problem that the classical 
phase space structure is typically unknown for a com- 
plex scatterer. A viable measurement protocol which is 
based solely on experimentally accessible quantities like 
the scattering matrix [l[ would thus be highly desirable. 

To resolve the contributions of individual phase space 
bands we present an approach based on the observation 
that all trajectories in the same band have a character- 
istic and very similar cavity dwell time. In analogy to 
eikonal theory we may thus "label" contributions from 
different bands by their respective dwell times (or path 
lengths) . In wave scattering the closest analogues to clas- 
sical dwell times are the "proper delay times", i.e., the 
eigenvalues of the Wigner-Smith time-delay matrix fl3j ] , 



Q = in——b = in 
oh 



rt r + ih rH' + iV 
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where the dots stand for the energy derivative 8e ■ Using 
the eigenvalues qi of Q to lift the unwanted degeneracy 



in the noiseless subspace is, however, non-trivial since Q 
has a different dimension (2N x 2N) than the transmis- 
sion matrix vt (N x N). Accordingly, the eigenstates 
of Q, in general, are scattering states injected from both 
leads, whereas the eigenstates of vt are injected from the 
left lead alone. As shown below, this mismatch is conve- 
niently resolved in the noiseless subspace where a basis 
of common eigenstates to both Q and tH can be found. 

Due to the hermiticity of Q its eigenstates \qi) form an 
orthogonal and complete set of states, to each of which a 
real "proper delay time" q t can be assigned. In the cor- 
responding matrix representation of this eigenproblem, 
= q i <f/ n , the 2iV-dimensional time-delay eigenvec- 
tors q™ = (q^l, (lifj) contain the complex coefficients of 
the eigenstates \qi) in the flux- normalized basis of incom- 
ing modes in the left (\n)) and in the right lead (|^')) : 
{%%) n = ( n h) and {$i%) n > = ("'I*)- Correspond- 
ingly, the outgoing coefficient vectors q° ut = (q"^, q"^ 1 ) 
contain the coefficients in the basis of outgoing modes: 
(%L)n = CH<?<> and (q™\, = (Tn'\ gi ), where T is 
the time-reversal operator of complex conjugation (T 2 = 
1 for spinless scattering). With q° ut = S q™ and S = S T 
for systems with time-reversal symmetry we define an 
anti-unitarity operator S = TS = S^T which maps the 
incoming coefficients of a time-delay eigenstate onto the 
incoming coefficients of the corresponding time-reversed 
state, (<fj OUt )* = Scfj 111 . As this operator S commutes 
with the time-delay operator, [S, Q]=0 (see [14(A), any 
non-degenerate time-delay eigenstate is time-reversal in- 
variant (up to a global phase e m , aeR): 5q* ln = e %a q^ n . 
For non-degenerate time-delay eigenstates whose incom- 
ing flux from one lead exits through both of the leads 
this time-reversal invariance implies that these states 
must also have incoming flux contributions from both 
leads. Such 2A^-dimensional time-delay eigenvectors q™ 
can thus not be reduced to an A-dimensional vector with 
incoming flux from the left lead alone. 

This restriction is lifted in the noiseless subspace, 
where the incoming flux from one lead also exits through 
just one of the leads. Consider a noiseless time-delay 
eigenstate with fully transmitted incoming flux from the 
left lead, tHq^ = q^, but no incoming flux from the 
right lead, qfy = 0. For this state, q™ = [^£,0], the 
commutator [H, Q] — implies that the time-reversed 
state, Sifj 111 = [0, {tq/1)*], is also a time-delay eigenstate 
with the same eigenvalue qi as g^ 1 ". Being a noiseless 
eigenstate of t'H' with incoming flux only from the right 
lead, H^ ln is clearly orthogonal to q™. We thus find 
that such "NOiseless Time-delay Eigenstates" (NOTEs) 
come in pairs of two which together form the basis of 
a doubly degenerate subspace associated with the time- 
delay eigenvalue qi. We emphasize that, in contrast to 
the common eigenbasis of 5 and Q in this subspace, 
[^£,±(^2)*] />/2, NOTEs are not time-reversal in- 
variant. Rather, two NOTEs forming a degenerate pair 
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are the time-reversed of each other like a classical trajec- 
tory and its time-reversed partner. Focussing now only 
on NOTEs injected from the left lead, we can determine 
their expansion coefficients q^l with Q from Eq. @, 
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For the last equality to hold, the following two condi- 
tions need to be fulfilled: (i) Q n q^f = q i g*™ and (ii) 

Q 21 q^l = 0. This central result of our Letter implies the 
following operational procedure to determine the expan- 
sion coefficients of NOTEs: In a first step (i) the eigen- 
states of the hermitian matrix Qn of dimension N x N 
are calculated. Out of this orthogonal and complete set 
of vectors the subset which, according to (ii) , lies in the 
null-space (kernel) of Q21, constitutes the desired set of 
common eigenstates of Q and vt. In practice, condition 
(ii) can be conveniently verified by a null-space norm 
Xi = IIQ21 QitW which determines the degree to which 
the normalized vector q^\ lies in the null-space of Q21. 
The quality of a NOTE should be the better the closer 
this measure x G [0, 00] is to zero. As the limiting value 
X —> is only reached for exact NOTEs with wavelength 
A —> 0, we need to test whether our approach works also 
for the realistic situation where A has a finite value. 

Consider, as a starting point for such a test, our pre- 
vious argument that NOTEs can only exist in the noise- 
less subspace with r = 0,1. We emphasize that this 
requirement, which can similarly not be fulfilled exactly 
for any finite value of A, does not explicitly enter condi- 
tions (i),(ii) from above. A good indicator for the validity 
of our approach is thus the degree to which NOTEs with 
finite values of A are, indeed, noiseless. For this purpose 
we calculate the eigenvectors q>\ of Qn [see condition 
(i)] and verify how the transmission of these states cor- 
relates with the corresponding null-space norm \i [ see 
condition (ii)]. We find that all eigenstates of Qn which 
closely fulfill the NOTEs condition (ii) of low %-values 
are, indeed, either almost fully transmitted, r w 1, or 
fully reflected, r « 0. This behavior is entirely absent 
in Random Matrix Theory (RMT) where no phase space 
bands exist (a comparison of the data for the rectangu- 
lar cavity with RMT is shown in Fig. [5^ and in [14|B 
for other scattering geometries). After these consistency 
checks we test whether NOTEs with very low null-space 
norms (x J$ 30, see Fig. [2^,) display in their wave func- 
tions the anticipated pronounced enhancements around 
individual bundles of classical trajectories. Our results 
based on the same scattering matrix data as for Fig. [TJd 
confirm that states with such low null-space norms x ai l 
feature highly collimated beam-like wave functions (see 
Fig.Eb-d and [l3]E). Quite different from arbitrary noise- 
less states (Fig. [T^,b), NOTEs feature Husimi distribu- 
tions that do not mix contributions from different phase 
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FIG. 2: (Color online) (a) Transmission T vs. null-space norm 
X for eigenstates of the matrix Qn in a rectangular cavity with 
different lead orientations (green dots). NOTEs with \ 
are noiseless and strongly deviate from RMT (red diamonds, 
see [l4|]B for details), (b)-(d) Wave function densities for 
NOTEs calculated with the same scattering matrix data as 
used for Fig. [TJd. As demonstrated by the Husimi plots in 
the bottom panels, each state is located on a single classical 
phase space band. Null-space projections \ are 00 4.7, (c) 
6.3, (d) 6.9. The insets in (c) illustrate the possibility to use 
NOTEs for transferring information between a sender (A) and 
a receiver (B) which by-passes a potential eavesdropper (E). 



space bands, thereby corroborating the successful opera- 
tion of our procedure. Without exception we find that in 
cases where NOTEs seem to feature contributions from 
more than one band (as in Fig.[2}i) all these bands belong 
to a single connected band in an extended zone scheme 
(like T1/T2/R2 in Fig. [It) jH. 

Our numerical results indicate furthermore that proper 
delay times of NOTEs do not only lift the degeneracy of 
noiseless states located on different bands, but that also 
the small dwell-time differences between trajectories of 
the same band do get increasingly well resolved in the 
limit A — > 0. Correspondingly, we find that the proper 
delay times qi of NOTEs on the same band are charac- 
teristically different from each other (rather than degen- 
erate). NOTEs thus fill individual phase space bands in 
a well-controlled fashion. Consider, e.g., the band Tl: 
starting from the state in Fig. [2]: the proper delay times 
and the transverse quantization of states on this band 
increase (see, e.g., Fig. [2jl) until, when the band is filled, 
the null space norm of states increases substantially (see 
[3|E), indicating a substantial overlap with phase space 
outside of the band. Such an increase in %-values is of- 
ten found to be accompanied by signatures of diffractive 
scattering at the sharp lead mouths (see [5jE). 

Since the operational procedure presented here does 
not rely on any specific assumptions concerning the type 
of scattering in a given system, we also applied it to more 
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FIG. 3: (Color online) NOTEs at k = 75.5n/d in cavities 
(a),(b) with a Sinai- type boundary shape and (c),(d) with a 
bulk disorder of amplitude Vo/E — 0.1 (see illustrations). The 
correlation of the disorder is long-range for (c) with correla- 
tion length kr c — 30n and short-range for (d) where kr c = 5n. 
Null-space norms x are (a) 1.7, (b) 24.8, (c) 4.9, (d) 65.7. 



complex scattering geometries. Consider first the Sinai- 
type billiard structure in Fig. [3^,b which features chaotic 
classical dynamics due to scattering at the circular part 
of the hard wall potential. We find that NOTEs which do 
not have any overlap with the curved part of the bound- 
ary (see Fig. have wave functions and proper delay 
times as in the rectangular cavity. For comparison, we 
also show in Fig.[3}D a state which bounces off the circular 
boundary. As reflected in its increased null-space norm 
(x = 24.8) the high instability of this state's bouncing 
pattern makes it much harder to resolve the (small) area 
of its phase space band. Consider next rectangular cavi- 
ties containing a static disorder potential with correlation 
length r c and an average amplitude V$. By studying such 
disordered cavities we find that our approach is restricted 
to the limit of weak and long-range disorder, Vq/E <C 1 
and kr c 3> 1 (in agreement with the validity criterion for 
the WKB / eikonal approximation [l8[ and previous work 
[Io|). Correspondingly, for the same fixed disorder am- 
plitude (Vo/E = 0.1) the z-shaped scattering state from 
Fig. [2)3 survives in the presence of long-range disorder 
(see Fig. [3}:), whereas no such state exists for short-range 
correlations (see Fig. |3Ji for the state with the closest 
time-delay value and [14jB,E for more details). 

We have performed additional tests (see 14(C) to ver- 
ify that all NOTEs which we find within the above 
limits are associated with individual classical trajec- 
tory bundles. This feature allows for a WKB/eikonal- 
type ansatz for the transmission matrix of NOTEs, t ~ 
Qi^R 1 e lSh ( E ^ h , in which the only part with a signif- 
icant energy-dependence is the action phase, Sb(E) = 
J b kdl, accumulated along bundle b. We show in [l4j]D 
that this ansatz fulfills the defining conditions for NOTEs 
(i),(ii) from above. 

We believe that our results open up many interesting 
possibilities for the experiment, where the cavities con- 
sidered here could, e.g., be an acoustic resonator (like a 



room) or an electromagnetic scatterer (like closely spaced 
buildings). In both these cases the collimated wave func- 
tions associated with NOTEs might allow to transfer 
information between a sender (A) and a receiver (B) 
such that the power to generate the signal is minimized 
and the transmitted signal is kept out of reach of an 
eavesdropper (E) (see illustration in Fig. [2J;). In this 
sense NOTEs offer clear advantages over arbitrary noise- 
less scattering states that do not display such beam-like 
wave functions in general (see, e.g., Fig. |TJj) . NOTEs 
may also have interestin g co nnections to phenomena in 
closed/decaying systems [15l4l7j |. 

In summary, we present an operational procedure 
for constructing scattering states which follow classi- 
cal bouncing patterns in coherent transport through 
cavities or complex scattering landscapes. In anal- 
ogy to WKB/eikonal theory we find that such ray- 
optical/beam-like scattering states arc determined by the 
condition of a fixed scattering time-delay. Our procedure 
is generally applicable to different types of wave scatter- 
ing (acoustic, electro-magnetic, quantum etc.) and relies 
solely on the knowledge of the scattering matrix. 

We wish to thank F. Aigner, J. Burgdorfer, A. Cresti, 
and A. Foerster for helpful discussions. Support by the 
WWTF and computational resources by the Vienna Sci- 
entific Cluster (VSC) are gratefully acknowledged. 
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A. TIME-REVERSAL OPERATOR 

In the main article we introduce the anti-unitary time- 
reversal operator S = TS which maps the incoming co- 
efficients of a scattering state onto the incoming coef- 
ficients of the corresponding time-reversed state. For 
time-reversal invariant systems like the ones we study 
the scattering matrix S is unitary symmetric such that 
the time-reversal operator S can also be written as fol- 
lows: S = S^T. This identity can be conveniently used 
to show that 5 commutes with the Wigner-Smith time- 
delay operator Q which, due to its Hermiticity, satisfies 
the following identities: Q = ih^ES^S — —ihS^dsS). 
In particular, we have: 

[Q,E] = {ih(d E S f )S} S f T- S f T{-ihS f (d E S)} = 
= ih{(d E S*)SS*T - S*ST(d E S)} = 
= ih{(d E S^)T - (d E S r )T} = . (E.l) 

B. STATISTICAL ANALYSIS OF NOTES 

Since we expect the presence of NOTEs in a given scat- 
tering system to depend on non-universal, system specific 
scattering processes, we should find no NOTEs in a sys- 
tem described by an entirely random scattering matrix 
as in Random Matrix Theory (RMT). Whereas RMT is 
thus not a suitable framework for studying NOTEs, the 
differences that appear between an RMT description and 
a full numerical solution of a non- universal scattering sys- 
tem will be very instructive for characterizing NOTEs in 
detail. 



Random Matrix Model 

For a suitable RMT model of scattering which pro- 
vides also the Wigner-Smith time delay matrix we follow 
the so-called "Heidelberg approach" reviewed, e.g., in [lj. 
This approach links the scattering matrix S of a system 
with the Hamiltonian H of the corresponding closed cav- 
ity through the coupling matrix W between the discrete 
bound states in the cavity and the continuum in the at- 
tached leads, 

S(E) = I-2^ E _ H l mW ^ W. (E.2) 



For time-reversal invariant systems the Hamiltonian ma- 
trix H is taken from the Gaussian orthogonal ensemble 
(/3 = 1). The mean spacing of the energy eigenvalues 
of H is adjusted to match the corresponding mean level 
spacing AE = 2ir/A of our cavities with size A. (Note 
that the chaotic cavities we study have a smaller scatter- 
ing area A and are thus compared with a different RMT 
ensemble than the regular and the disordered cavities.) 
The mean absolute value of the random coupling matrix 
elements of W is adjusted such as to produce a scattering 
matrix S with, on average, balanced total transmission 
and reflection, (T) sa (R) sa (T') « (R') ^ Af/2 |. 

For calculating the Wigner-Smith time-delay matrix 
based on Eq. (|E.2I) . Q = ih(d E S^)S, we make the usual 
assumption |l| that the coupling matrix elements in W 
are only weakly energy-dependent such that the deriva- 
tive dsS^ affects only the explicit energy dependence in 
the denominator of Eq. (|E.2|) . To check our results we 
have verified that the Wigner-Smith time-delay matrix 
Q obtained in this way gives rise to eigenvalues qi which 
are statistically distributed according to the RMT pre- 
dictions put forward in [3] (not shown). 

Comparison of numerical data with RMT 

It is now of interest to compare the results obtained 
with this RMT approach with those from our numerical 
calculations for the regular, chaotic, and disordered cavi- 
ties. The RMT data was produced based on an ensemble 
with altogether 20 random scattering matrices following 
Eq. (|E.2j) . To obtain a comparatively large ensemble of 
statistical data also for the cavities, we performed the 
numerical wave calculations with 20 different disorder re- 
alizations (for the disordered cavities) and with 20 differ- 
ent lead configuations (for the regular and the chaotic 
cavity). In the latter case, we shifted the right lead in 
equidistant steps from the bottom to the top position 
(the left lead was kept fixed at the top position). To 
counterbalance the increased numerical effort necessary 
for calculating such ensembles of cavities, we reduced the 
wavenumber from k = 75.5tt/cI (as used for the plots in 
the main text) to k = 45.5w/d. 

In both the RMT and the cavity calculations we fo- 
cused, in particular, on the eigenstates q" L of Qn which, 
according to condition (i) in our procedure (see main 
text), are the key ingredients for NOTEs. In Fig. E.I we 
plot the total transmission of these states as a function of 
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Fig. E.l: Transmission versus null-space norm \ f° r eigen- 
states of the submatrix Qn of the Wigner-Smith time-delay 
matrix Q (wave number k = 45.5n/d). (a) Schematic of the 
results shown in Fig. 2a (main text) to illustrate the crossover 
between particle/classical and wave/quantum scattering, (b)- 
(d) Same as Fig. 2a (main text) for (b) the Sinai-type cavity 
with chaotic classical dynamics as well as for weakly disor- 
dered cavities (Vb = 0.1E) with (c) long-range (kr c — 15n) 
and (d) short-range (kr c = 2.57r) disorder. The red diamonds 
show the data obtained from the RMT calculations from alto- 
gether 20 randomly assembled S-matrices (for better visibility 
we use an increased symbol size for the RMT data). Strong 
deviations between the cavity and the RMT data are observed 
which point to the presence of NOTEs with small null-space 
norms \ ~ ► and noiseless transmission T(l — T) — > 0. The 
presence of a disorder potential (in particular, with short- 
range correlation) shifts the cavity data towards the RMT 
distribution and strongly reduces the amount and the quality 
of NOTEs. 



the corresponding null-space norm x — ||Q2iSj 1 £,|| which, 
according to condition (ii) measures whether a state can 
be identified as a NOTE. In this plot already a number of 
important features become apparent: the data from the 
RMT calculations (red diamonds) accumulate in a very 
restricted region of the parameter space: the norm \ of 
the RMT states stays within a restricted interval of x and 
their transmission clusters very strongly around the mean 
transmission T — 0.5. This situation is quite different for 
the cavities: whereas the upper bound of the norm \ is 
comparable to the corresponding RMT value, the lower 
bound of x is more than two orders of magnitude smaller 
than in RMT, suggesting a strongly non-universal be- 
havior. Furthermore we observe that states with very 
low values of the norm x are a U clustered around the 
noiseless transmission values very close to T = and 
T = 1. 

This observation demonstrates that, indeed, eigen- 
states of Qn with a vanishing nullspace norm [x — > 0] 
have noiseless transmission [T(l — T) — > 0]. Equivalently, 
each scattering state injected from just the left lead [i.e., 
with Qiji = 0] which is a time-delay eigenstate is found to 
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Fig. E. 2: Null-space norm \ — WQiiQiXW versus noiseless 
norm a = 4T(1 — T) for the four different cavity types in- 
vestigated. The same data set and color coding is used as in 
Fig. E.l. A strong positive correlation between the two norms 
X and a is visible for all data points that are strongly non- 
universal, i.e., for parameter values that are different from 
those where the RMT results cluster (see red diamonds). 



be either fully transmitted or fully reflected. NOTEs are 
thus, indeed, noiseless. We can further write the trans- 
mission T of a NOTE as T = J2i \ c i\ 2 T~i, where the r, 
are the transmission eigenvalues of Pt and the Ci are the 
expansion coefficients of the NOTE in the eigenbasis of 
tH. From this we conclude that the noiseless transmis- 
sion values of NOTEs, T « (T w 1), imply that NOTEs 
can only have non-negligible expansion coefficients Cj for 
those transmission eigenvalues Tj which are themselves 
noiseless, Tj s» (rj ps 1). NOTEs can thus only be 
composed of noiseless transmission eigenchannels. 

We emphasize that the strong correlation between 
time-delay eigenstates and noiseless scattering (which 
does not explicitly enter our construction of NOTEs) per- 
sists here in the limit of finite channel numbers where 
neither the null-space norm x is exactly nor the trans- 
mission T is exactly noiseless for any state. To highlight 
this strong correlation we show in Fig. E.2 the null-space 
norm x of Qu eigenstates as a function of a = 4T(1 — T), 
which measure for noiseless scattering can vary between 
and 1. For all states which display non- universal be- 
havior (deviating from the RMT-result) a strong posi- 
tive correlation between the two measures x an d f is 
observed, suggesting that we can use these two measures 
interchangeably for testing the quality of a NOTE. We 
have checked explicitly that the statistical spread in the 
correlation between x an d cr (see Fig. E.2) is further re- 
duced when considering states on the same phase space 
band. 

The statistical data shown in Fig. E.l and Fig. E.2 
also demonstrate that NOTEs are completely absent in 
the RMT description along Eq. (|E.2[) . This observation 
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Fig. E.3: Eigenvalues qi of Qn versus noiseless norm a of the 
corresponding eigenvectors q/1 . The same data set and color 
coding is used as in Fig. E.l. As expected from our model, all 
states with a small noiseless norm a correspond to short delay- 
time values qi. A clustering of these qi -values around the 
classical values belonging to individual phase space bands is 
clearly observed. Note, e.g., that a delay time of 100 (200) cor- 
responds to a single (double) traversal of the cavity along the 
upper cavity boundary in those cavity configurations where 
the right lead is (is not) located at the top-most position. The 
maximum delay-time to where such non-universal deviations 
with a < 0.1 occur are indicated by the black horizontal bars. 
The values of to are reduced in chaotic vs. regular ballistic 
cavities and in short-range vs. long-range disordered cavities. 



is perfectly in line with the expectation that the presence 
of NOTEs requires extended bands in phase space. As 
such non-universal features are clearly absent in the RMT 
description, this model does not give rise to any NOTEs. 
Rather, the eigenstates of the random matrix Qn do not 
constitute a preferred basis in terms of transmission (as 
opposed to the eigenstates of t't). Their transmission 
values thus cluster around the average value T ss 0.5 and 
are strongly suppressed near the noiseless values T w 
0, 1. As discussed in the main part of the paper, one 
can show by employing the time-reversal operator 3 that 
such states which lie outside the noiseless subspace can 
never be time-delay eigenstates when injected from just 
one of the leads. Correspondingly, the Q21 null space 
norm \ of these RMT states always remains larger than 
the (small) %-values of NOTEs. 

From the above results we conclude that the noiseless 
eigenstates of Qn are NOTEs for which the correspond- 
ing eigenvalues qi of Qn are time-delay values. Following 
our model we expect that these time-delay eigenvalues 
correspond to the time spent in the cavity when scat- 
tering along a given trajectory bundle (see section C for 
more details). Consequently, all states that belong to the 
same bundle should feature similar values of qi, which we 
expect to lead to a clustering of ^-values around the clas- 
sical time delay values of those trajectory bundles which 



can be fully resolved by a scattering state. As, in turn, 
the degree of this resolution is measured by the noise- 
less norm a (or, alternatively, by the null-space norm 
x) it is instructive to plot the eigenvalues qi of Qn as 
a function of these norms a (or \). The resulting plot 
shown for the noiseless norm a in Fig. E.3 nicely con- 
firms the expected behavior: we first observe that no 
states with both very long time-delays and small a- values 
are found, in correspondence with the expectation that 
NOTEs can only form on short-lived/stable phase space 
bands. States with small norm values a show the most 
pronounced deviations from RMT and correspond exclu- 
sively to short-lived states with time-delay values on the 
lower end of the entire distribution. For a single cavity 
these non-universal (^-values cluster along horizontally 
elongated regions, corresponding to states on the same 
phase space band. For such states the time-delay val- 
ues qi are similar, but the norms a vary considerably, 
depending on the overlap of a state with regions outside 
the band. In Fig. E.3 such a clustering of values is still 
visible but blurred due to the different lead configura- 
tions that enter the statistical data set. Those NOTEs 
out of this set with the lowest values of x correspond to 
scattering states which propagate solely along the upper 
cavity boundary where the focussing is enhanced by the 
top cavity wall. In contrast, all of these non-universal 
features are drastically reduced for states with a norm a 
near the maximum value 1. 



Signatures of regular /chaotic/disorder scattering 

Since NOTEs are highly non-universal and system- 
specific scattering states, their number and their qual- 
ity should depend significantly on the type of scattering 
present in a given scattering region. To investigate this 
issue in more detail we calculated all of the above plots 
for the regular, chaotic and disordered cavities, respec- 
tively. As a measure for the influence of these different 
scattering mechanisms we introduce a critical time scale 
To which represents the maximal time up to which time- 
delay eigenvalues of a given quality (do = 0.1) exist for a 
given system. As this time scale To measures for how long 
scattering states follow the classical bouncing patterns of 
NOTEs (see section C) t can loosely be associated with 
the Ehrenfest time (see Refs. [5-10] in the main paper). 

In all of the above plots we consistently find that the 
strongest deviations from RMT occur for the regular cav- 
ity which, correspondingly, also features the largest value 
of To (see black horizontal bars in Fig. E.3). Note that 
for the regular cavity deviations from RMT do persist to 
a certain degree also for states with time-delays longer 
than T , indicating that in regular cavities also very long- 
lived states may carry non-universal features. Compar- 
ing these results with those for the chaotic cavity we ob- 
serve the following interesting behavior: whereas states 
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with the shortest time-delays feature very similar values 
of the norm a in the regular and the chaotic cavities, 
t is considerably shorter for the chaotic scatterer than 
for the regular one. This can be intuitively explained 
by the fact that short-lived NOTEs which do not have 
any overlap with the circular part of the boundary in the 
chaotic dot are the same as the corresponding states in 
the regular cavity. Those very stable states correspond- 
ing to the smallest a values thus appear in both cavity 
types with similar parameter values. Longer-lived states 
do eventually also explore the circular boundary part in 
the chaotic cavity which leads to increased wave spread- 
ing and thus to a reduced value of To . In the cavity with 
long-range disorder the situation is again conspicuously 
different from the chaotic cavity, although both types of 
cavities give rise to classical chaotic scattering: in the 
long-range disordered cavity even the most short-lived 
states are drastically affected by the disorder since it is 
contained in the bulk and affects all scattering states (in 
contrast to the chaotic cavity where the circular bound- 
ary part can be avoided) . Correspondingly, the minimum 
value of the norm a is much higher for the disordered 
cavity than for both the regular and the chaotic cav- 
ity. When changing the correlation length of the disor- 
der to smaller length scales (with the disorder amplitude 
remaining unchanged at Vq/E = 0.1) the effect of the dis- 
order is further increased: here even the states with the 
shortest delay time have a values which are three orders 
of magnitude larger than in the ballistic cavities without 
disorder. Also To is here much smaller than in all other 
cavity types. This behavior can be attributed to the ef- 
fect of stochastic scattering which is induced whenever 
the correlation length of the disorder is of comparable 
magnitude or smaller than the wavelength. 

NOTEs and general time-delay eigenstates 

In the main part of the article we argue that NOTEs 
should appear as a special subset of eigenstates of the 
Wigner-Smith time-delay operator Q. A convenient way 
to test this relation between NOTEs and general time- 
delay eigenstates is to compare the eigenvalues and eigen- 
vectors of the submatrix Qu (which we use to gener- 
ate NOTEs) with those of the general Wigner-Smith 
time-delay matrix Q. In particular, we expect that the 
eigenvectors and eigenvalues belonging to a NOTE (i.e., 
for which x,a — ¥ 0) should emerge for both matrices, 
whereas the remaining eigenstates (with higher values of 
X, cr) may again be close to the corresponding RMT val- 
ues for each of these two matrices. To check whether this 
expected behavior is, indeed, fulfilled consider Fig. E.4: 
this plot shows the eigenvalues qi as a function of the 
noiseless norm a obtained here from the entire matrix Q 
rather than only from Qu as in Fig. E.3. The two sets 
of eigenstates from Q and Qu are of different size since 
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Fig. E.4: Same as Fig. E.3 but with the time-delay eigenval- 
ues qi and the noiseless norm a determined from the eigen- 
problem for the whole Wigner-Smith matrix Q (see text for 
details). As predicted by our model, a clear correspondence 
with Fig. E.3 can be observed for small values of a. This indi- 
cates that NOTEs as identified by the procedure proposed in 
the main text naturally emerge as eigenstates of the Wigner- 
Smith time-delay matrix. 

the eigenvectors q™ of Q contain the incoming mode 
coefficients from both leads, q™ = (q"l, Qi%), not just 
from the left lead alone. We thus calculate the norm 
a — 4T(1 — T) for Fig. E.4 just based on the transmis- 
sion T associated with the (normalized) injection coeffi- 
cients from the left lead, i.e., with qi'l/Wq^lW- Since we 
expect NOTEs to appear in the eigenvectors of Q as a 
near degenerate doublet (see main text) the normaliza- 
tion removes here the effect of mixing between these two 
degenerate states. Figure E.4 obtained in this way dis- 
plays the same characteristic and non-universal features 
as Fig. E.3. The way in which these features are similar 
in both plots corroborates our approach. 

The small differences between Figs. E.3, E.4 which re- 
main for small values of a can be well described with the 
help of degenerate perturbation theory (not shown). In 
this description the "perturbation" is given in terms of 
the off-diagonal matrices Q12, Q21 which mix the eigen- 
states of the "unperturbed" matrices Qu, Q22- Eigen- 
states and eigenvalues of Qn thus get perturbed through 
admixtures from other states which, in further conse- 
quence, lead to the above differences. In turn, these ad- 
mixtures should vanish when the perturbation of a Qn- 
eigenstate q^l by the matrices Q12 or, equivalently, by 
Q21 goes to zero. Note that this situation is exactly real- 
ized when the null-space norm x = IIQ219/1JI (or, equiva- 
lently, the noiseless norm a, see discussion above) goes to 
zero. This reasoning not only substantiates the fact that 
the norms x an d are a suitable means to measure the 
quality of a NOTE but also that in the numerically un- 
reachable limit of infinitely many lead modes (N — > 00) 
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we should expect the Wigner-Smith time-delay matrix to 
feature exactly degenerate pairs of NOTEs. 

C. NOTES AND CLASSICAL TRAJECTORY 
BUNDLES 

An interesting question to ask is whether all NOTEs we 
find are associated with classical trajectory bundles and 
vice versa. To answer this question we have performed a 
number of tests: our work builds on the established re- 
sult from the literature (see, e.g., |4|) according to which 
phase space bands with a size larger than Planck's con- 
stant give rise to noiseless scattering channels. In the 
weakly disordered regime (which we consider) each of 
these bands can be associated with a bundle of classi- 
cal trajectories of equivalent topology and similar length 
(a situation which is in general not fulfilled for strongly 
disordered systems). To verify whether all NOTEs are 
associated with such classical trajectory bundles we thus 
need to test if the NOTEs generated with our procedure 
[see conditions (i),(ii) on page 3 in the main paper text] 
are composed exclusively of noiseless scattering channels. 
As demonstrated in section B (see, in particular Figs. E.l, 
E.2) all NOTEs which we find in the ballistic structures 
(regular and chaotic), indeed, have this property. We also 
find that in the weakly disordered structures where such 
trajectory bundles are disturbed, NOTEs are of consider- 
ably lower quality. Furthermore, the quality and quantity 
of NOTEs correlates strongly with the area of a given 
phase space band: consider, e.g., that no NOTEs are 
found with very large time-delays (see Figs. E.3, E.4) as 
these are associated with very small phase space bands. 
We conclude from these observations that all NOTEs we 
find can, indeed, be associated with classical trajectory 
bundles. 

The above arguments still leave open whether each 
NOTE can be associated with only a single or multiple 
trajectory bundles. Since NOTEs are time-delay eigen- 
states and different bundles have, in general, different 
time-delays our procedure should yield only a single bun- 
dle per NOTE already by construction. To test whether 
this expectation is, indeed, fulfilled we performed the 
following additional tests: We have checked that each 
time-delay eigenvalue qi associated with a NOTE in the 
regular cavity corresponds to a classical time-delay value 
from a single classical phase space band: Based on the 
classical dynamics in a rectangular cavity the product 
of a proper delay time qi with the eigenstate's average 
longitudinal velocity, (\v l x \) — h(\k x \)/m, must be an 
integer multiple of the cavity width W , corresponding 
to the number M of the state's left-right cavity traver- 
sal, qi |u* | = MW. With the average longitudinal mo- 
mentum (\k x \) of state i being determined as follows, 

(141) = E,ti l(0)JVfc 2 -(Wtf) 2 , we find this cri- 
terion to be fulfilled with a relative error below 2% by all 



the NOTEs for which \ 20. This result confirms that 
for the rectangular cavity our numerical procedure yields 
proper delay times in very good agreement with classical 
expectation. In the chaotic cavity we have checked that 
NOTEs which have no overlap with the circular bound- 
ary feature a one-to-one correspondence with the equiv- 
alent states in the regular cavity (both in terms of the 
time-delay eigenvalue and in terms of the wave function). 
For the remaining NOTEs in the chaotic cavity (which 
do have an overlap with the circular boundary) and for 
all NOTEs in the disordered cavities we have checked 
explicitly with the help of the corresponding scattering 
wave functions that these NOTEs all feature highly col- 
limated wave function patterns as expected for states 
that are supported by classical trajectory bundles. De- 
tails can be found in section E, where we have provided 
the wave function patterns of NOTEs on specific phase 
space bands (in the ballistic cavities) and all NOTEs with 
X < 30 in the disordered cavities. 

Whereas the above tests suggest that each NOTE 
we found is associated with a single phase space 
band/trajectory bundle this association does not always 
hold in the opposite direction: as stated above, NOTEs 
can only be formed on phase space bands with a size 
larger than h (a condition which we have explicitly ver- 
ified in an earlier paper Q). A further restriction to 
form NOTEs on such noiseless channels comes from the 
well-defined time-delay associated with NOTEs. This re- 
quirement imposes an additional constraint which makes 
it harder for NOTEs to fit on a given phase space band 
as compared to the transmission eigenstates of t't. Cor- 
respondingly, we find that the transmission of NOTEs 
is generally not as noiseless as that of the transmission 
eigenstates (not shown). 

We mention, parenthetically, that imperfections both 
in the formation of noiseless channels and NOTEs are 
mostly reflected by those states labeled as "intermedi- 
ate" in Fig. E.la: these states deviate from RMT statis- 
tics but do not have a good quality as NOTEs. Such 
intermediate states are found to persist in the presence 
of weak disorder (Fig. E.lc,d). 

D. WKB/EIKONAL ANSATZ FOR NOTES 

In the main part of the paper we showed that NOTEs 
are defined by the two conditions (i) Qn q^\ = qi <?*™ an d 
(ii) Q21 <f/£ = 0. In this part of the appendix we inves- 
tigate how these equations can be satisfied, based on the 
observation from section C that NOTEs are time-delay 
eigenstates located on individual phase space bands. We 
assume, for simplicity, that we are dealing with a trans- 
mission band, leading to a fully transmitted NOTE. Our 
arguments can, however, be applied to reflection bands 
as well. 

Consider first the eigenproblem in (i): Qn q"l = 
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ih(f^r + tH) q^l — qiq"l- Using the property 
that q^l is a fully transmitted noiseless state, 
t tQiL — QiL> we have r^rq^l — f* = 0, leaving 
us with Qnq^l = ihtHq^l — qiq^l for the evaluation 
of which we perform a singular value decomposition 
(SVD) of t = with ti = UiO-ivj. Here the 

right-singular vectors Hi = r ; ln are the eigenvectors 
of tH, the left-singular vectors Ui = rrV ln ) are the 
eigenvectors of ttf = (t'H'}*, and the singular values 
are the square roots of the corresponding transmission 
eigenvalues, cr, = ^ffi € [0, 1]. In wave transport through 
complex scattering landscapes the non-deterministic 
singular values < Ui < 1 and their corresponding 
singular vectors v\ and u, feature a very strong energy 
dependence as reflected, e.g., in the familiar conductance 
fluctuations. In the degenerate noiseless subspace where 
Ui = 0, 1 the situation is more complicated: states which 
are noiseless due to an appropriate interference of path 
contributions (like in Fig. QJi) also feature a significant 
energy dependence induced by multi-path interference. 
For those states that are noiseless due to a resolution 
of the underlying classical phase space (like in Fig. [T]d) 
the singular values <Tj all stay at the energy-independent 
degenerate values of <Tj =0,1 (see, e.g., [4j). Since the 
corresponding singular vectors are, typically, distributed 
over several transmission (reflection) bands for a{ = 1 
(<Tj = 0) their energy-dependence is not well defined 
unless we fix this distribution. As we argue above, a 
distribution where each state is located on a single band 
can be achieved by demanding that the noiseless vectors 
Vi are NOTEs, Vi — and Ui — tq"l. On a single 
band all involved trajectories will accumulate a very 
similar phase, leading to energy-independent construc- 
tive interference such that the incoming/right singular 
vectors have only a very weak energy dependence, 
BeVi ~ 0. In this relation the arbitrary global phase that 
we are free to choose in any SVD was fixed for the vector 
Vi such as to minimize its energy-dependence. This 
condition, in turn, also determines the corresponding 
outgoing/left singular vector, Ui = M,°e*' £ ', to consist 
of an energy-independent part, Oeu® s» 0, and an 
energy- dependent phase shift, e l ^ E \ The phase shift 
4>b(E) — Sb{E)/h— 7r^b/2 can be associated with the 
action Sb(E) = J b kdl and the Maslov index /ib accumu- 
lated along the phase space band b. With the transmis- 
sion amplitudes for the transmitted NOTEs now being 
ti = 9j I fl*e I<?i ' , g i 1 ^, we can further simplify condition (i): 

Qu0, = mhq™ = ih[d E (thq^ L )-^a E (t^ L )] = 



(dESb) q"l, from which we obtain the time-delay eigen- 
value Oj = d E Sb, in full analogy to the one-dimensional 
case [6(. Similarly, condition (ii) reads: Q21 ^'l = 
ihr'Hq™ = iK[d E {r'hq™)-r*d E {tq>X)] " = 
ih \pE{r'^t q^D — (dESb) r'Hq^] , which expression 
must be zero since the unitarity of the S'-matrix implies 
r'H = —t''r and t'^rq^l = for a noiseless transmission 
state. 

We have thus verified that WKB / eikonal-type scatter- 
ing states which are characterized by their phase shift 
accumulated along a given trajectory bundle satisfy the 
two defining equations (i),(ii) which we have put forward 
to construct NOTEs. 



E. WAVE FUNCTION IMAGES 

In the figures attached below we provide a number of 
wave function images as obtained for the four different 
cavity types discussed in the main article. In the regu- 
lar and chaotic cavity we show all states belonging to a 
specific phase space band: Figs. E.5, E.6 illustrate how 
NOTEs systematically increase their transverse quanti- 
zation until, when a band is filled, overlap with phase 
space regions outside this band lead to a significant drop 
in the quality of NOTEs. For the disordered cavities we 
display all wave function patterns with a quality below 
X < 30: Figs. E.7, E.8 show that all states of this qual- 
ity still feature very collimated wave functions in spite of 
the presence of weak disorder. Note that in the cavity 
with long-range correlations in the disorder the number 
of NOTEs of this quality is much higher than with short- 
range disorder. 
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Fig. E.5: Wave function densities of NOTEs in the rectangular cavity at wave number k — 75.5-7r/d: States associated with the 
transmission band T2 (see Fig. lc in the main text) are shown with their respective null-space and noiseless norms \y a - With 
increasing transverse quantization the quality of states (as measured by \, a) deteriorates monotonically. This decrease in the 
quality of NOTEs is accompanied by the appearance of diffraction effects at the sharp lead mouth (see arrow in the bottom 
row). 




Fig. E.6: Wave function densities of NOTEs in the chaotic cavity at wave number k = 75.57r/d: States associated with the direct 
transmission path are shown with their respective null-space and noiseless norms \i a - With increasing transverse quantization 
states have more overlap with the circular boundary part. For low-quality NOTEs diffraction effects are observed also here (see 
arrow in the bottom row). 
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Fig. E.7: Wave function densities of NOTEs in the disordered cavity with long-range disorder at wave number k — 75.57r/d: 
all NOTEs with a null-space norm x J; 30 are shown in the order of their time-delay value qi. The number of states with this 
quality is here much larger than in the case of short-range disorder with the same disorder amplitude, Vo — 10%E (see Fig. E.8 
for comparison). The employed potential landscape with disorder correlation length kr c = 307r is shown in the rightmost panel 
in the bottom row. 




Fig. E.8: Wave function densities of NOTEs in the disordered cavity with short-range disorder at wave number k — 75.5n/d: all 
NOTEs with a null-space norm % ^ 30 are shown in the order of their time-delay value qi. The disorder amplitude (Vo = 10%E) 
is here the same as in Fig. E.7, but the number of states with quality x 30 is strongly reduced as compared to the case of 
long-range disorder. The employed potential landscape with disorder correlation length kr c = 5tt is shown in the rightmost 
panel in the bottom row. 



